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Abstract 



We present a 2+1 dimensional quantum gauge theory model with correlated 
"^ ' fermions that is exactly solvable by bosonization. This model gives an effective de- 

scription of partially gapped fermions on a square lattice that have density-density 
interactions and are coupled to photons. We show that the photons in this model are 
massive due to gauge-invariant normal-ordering, similarly as in the Schwinger model. 
'^ ■ Moreover, the exact excitation spectrum of the model has two gapped and one gapless 

^ , mode. We also compute the magnetic field induced by an external current and show 

Q \ that there is a Meissner effect. We find that the transverse photons have significant 

effects on the low-energy properties of the model even if the fermion-photon coupling 
is small. 

^ ' The possible violation of Landau's Fermi liquid theory in models of strongly interacting 

Q ■ fermions has been an actively researched problem for many years. Interest in this topic 

quickly grew with the discovery of the cuprate high-temperature superconductors in 1986 

O \ pLj , and the realization that these materials display many properties not described by Fermi 

liquid theory [2j. Early on, it was suggested that models of Hubbard-type capture the 
strongly correlated physics of cuprates [21 HI El E]. While it has proven very difficult to 
do reliable computations for two dimensional (2D) such models, an excellent theoretical 
understanding has been obtained for the corresponding one-dimensional (ID) cases. This 

5^ \ understanding is largely based on models that are exactly solvable. A famous example is 

the Luttinger model [71 E], which provides an effective description of interacting fermions 
on a ID lattice, and which has become a prototype of non- Fermi- liquid behavior [9l [TO] . 

It has been known for quite some time that fermions coupled to dynamical photons 
can have non- Fermi- liquid behaviour [HI [12]. One argument against this mechanism being 
relevant for real materials is the smallness of the fine-structure constant (a ~ Y37)) which 
governs the strength of interactions between matter and transverse photons. However, 
various scenarios have been proposed and explored in which effective photon-like gauge 
fields arise in the low-energy limit of models for strongly correlated fermions [T3l [TH [15] ; 
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see also [10]. In these instances, the effective couphng constant need not be smalL The 
computations to explore this mechanism in 2D are usually based on approximations that 
are difficult to justify and, again, things are much better understood in ID due to the 
existence of exactly solvable prototype models. For example, the (1+1)D quantum gauge 
theory obtained by minimally coupling the Luttinger model to dynamical photons is exactly 
solvable [16] . This model is a generalization of the Schwinger model [17] and, similarly as in 
the latter model [18], the photons have a non-zero mass that arises due to gauge- invariant 
normal ordering. 

Despite a concerted effort, much remains to be understood when it comes to strongly 
correlated fermions in 2D. We believe that there is a need for models that are amenable to 
exact computations. One pioneer in this direction is Mattis who, already in 1987, presented 
a 2D interacting fermion model that can be mapped exactly to non- interacting bosons [IH] . 
Despite its potential, Mattis' model did not receive much attention (notable exceptions 
are [20| |2T| [22]). In recent work, we derived an effective model describing fermions on 
a square lattice with local interactions and close to half filling [231 [2H [25]. Since this 
model is similar to Mattis' [19], we refer to it as the Mattis model. In [26] we obtained a 
complete solution of the Mattis model, including all Green's functions. We also proved that 
the fermion two-point functions have algebraic decay with exponents that depend on the 
interaction strength, which is a hallmark of Luttinger liquids P]. The Mattis model is not 
relativistically invariant but, as already noted in [19], this is natural for effective models of 
square-lattice fermions. 

In this Letter, we present a (2-|-l)D quantum gauge theory obtained by coupling the 
Mattis model to a dynamical f/(l) gauge field. We refer to the latter as photon field, even 
though other physical interpretations are possible. Our main result is that this model 
is exactly solvable by bosonization. We also show that this model has some remarkable 
properties, including a non-zero photon mass due to gauge-invariant normal ordering, and a 
Meissner effect. We find that there is a renormalization of the bare charge, which determines 
the strength of fermion-photon interactions. Furthermore, the gauge theory contains a 
gapless mode that is strongly affected by the presence of the dynamical photons even if the 
fermion-photon coupling is weak. The idea of mapping fermions to free bosons has been 
used before to study 2D interacting fermions coupled to gauge fields [271 [2H], but no exactly 
solvable fermion model was clearly identified before. 

{Notation: /i, i/ = 0,+,— are space-time indices, s,s' = +, — space indices, r,r' = ± 
chirality indices; the space-time metric signature is (—,+,+); x = (ct,x) are space-time 
coordinates with 2D positions x = {x~^, x~); c is the velocity of light; dg = d/dx'^ are spatial 
derivatives; iprj are standard fermion field operators; A^ is the gauge potential. Eg are the 
electric field components, and B = d+A^ —d-A+ is the magnetic field. Common argument 
X of field operators are suppressed whenever possible.) 

The relation of the Mattis model to square-lattice fermions is described using Figure [H 
which shows the Brillouin zone corresponding to a square lattice (dashed large square) 
divided into regions of non-equal sizes. Close to half filling, mean field theory predicts that 
the system is partially gapped, and there is an underlying Fermi surface in the so-called 
noda/ regions (the four tilted rectangles) [26j. We model this underlying surface by straight 




Figure 1: The Brillouin zone corresponding to a 2D square lattice with lattice constant a. 
The quantum gauge theory provides a low-energy description of fermions with momenta 
in one of the four tilted rectangles labeled by (r, s), r, s = ±. Fermion degrees of freedom 
outside the rectangles are assumed to be gapped. 



arcs, which either corresponds to a truncated Fermi surface or the portion of a closed Fermi 
pocket having dominant momentum occupation. The Mattis model describes the fermion 
degrees of freedom in the vicinity of these arcs. It is written in terms of four fermion field 
operators il^r,si ''", -s = ±, in one-to-one correspondence with the nodal regions. The quantum 
field theory limit making this model amenable to bosonization amounts to removing the 
momentum cutoff orthogonal to the arcs (indicated by the arrow in the nodal (+, +)-region, 
for example), which is possible after normal ordering |24|. The arc-picture underlying our 
derivation of the Mattis model is supported by renormalization group studies of weakly 
coupled 2D Hubbard-like systems [29j . It is also a signature feature of the pseudogap phase 
as observed in angle- resolved photoemission experiments on hole-doped cuprates |30j . 



The Hamiltonian of the Mattis model is 






(1) 



with {V'rs(^)' V'r' s'(y)} ~ ^r,r'^s,s'^'^{^ — y) i ©tc, and colous denoting standard fermion 
normal ordering with respect to a Dirac vacuum [26] . These definitions are formal since an 
important UV regularization has been suppressed: while the Xs-component of x in iprj (x) is 
continuous, the x_s-component is discretized to integer multiples of a UV cutoff a, and thus 
the integrals and Dirac deltas have to be interpreted as partial Riemann sums and Kronecker 
deltas [26]. The coupling constants scale with the UV cutoff as gr,s,r',s' = anvF(^^i6s,s'Sr-r' + 
725s,_s'/2), with the Fermi velocity vp > and dimension-less constants 71^2 such that 
I71I < 1 and I72I < |1 + 7i|. The above scaling of the coupling constants is not only 
obtained by deriving the model from lattice fermions [2l|, but it also ensures that Hm has 
a non-trivial limit as a — )■ 0. The restrictions on 71 2 are to ensure stability of the Dirac 
vacuum 



The Hamiltonian of our quantum gauge theory model is obtained by couphng the Mattis 
Hamiltonian in ([1]) to a dynamical electromagnetic field, 



r,s=± r',s'=± /o") 

,2 d2 \ X 



with [yls(x), i?s'(y)] = i^s,s'<^^(x — y), etc., eo the bare charge, x ■ ■ ■ x boson normal ordering, 
and o ■ ■ ■ o a gauge-invariant generalization of fermion normal ordering (see below). The 
Gauss law operators generating gauge transformations A^ — )■ A^ + dsX^ etc., are 

G\x\ = [d'xxYl (-^^^« + ^oYl °^i.^r,. ^) , (3) 

"^ s=± r=± 

and the physical states are those annihilated by G[x], for arbitrary real- valued functions 
x(x). Note that, except for the normal-ordering procedures, the Hamiltonian in ([2]) is 
obtained from the Mattis Hamiltonian by standard minimal coupling: introduce the action 
corresponding to the Mattis Hamiltonian in ([T]); couple to Abelian gauge fields using the 
substitution —id^ — )■ —id^ + CqA^; add the Maxwell term c^F^^F'^^ /A] perform the usual 
Dirac procedure for systems with constraints [31] to obtain the formal Hamiltonian (i.e. 
without normal ordering) and the formal Gauss law constraint. We note in passing that 
the gauge field operators are well-defined without UV regularization (details will be spelled 
out elsewhere). 

As already indicated, fermion normal ordering plays a key role for our model. Indeed, 
formulating a sensible quantization of the classical theory is non-trivial, even in the absence 
of gauge fields: In order for the Hamiltonian to be bounded from below (have a ground 
state), we need to normal-order all fermion bilinears with respect to the Dirac vacuum 
in which all negative energy states are occupied. An important consequence of normal- 
ordering is that the fermion densities Jr^s = '■ i^l si^r s '■ obey the anomalous commutator 
relations [SSj 

[Jr,six), Jr',s'{y)] = r5ry5s,s'{'2'Kia)~^ds5'^{x - y), (4) 

and the non-interacting part of ([T]) can be expressed in terms of these densities using the 
operator identity 

d^X -.tpl gV {—ids) tpr s- = '^(^ / d^^ xJ^gX (5) 



(see Propositon 2.1 in [26] for proofs of these statements). The commutation relations in 
(jl]) imply that 

ds^s = VTTd{J+,s + J~,s), n^ = V7l5.{-J+^s + J~,s) (6) 

define boson operators obeying the usual canonical commutator relations [$s(x), Hs/(y)] = 
iSs^s'S'^i'^ — y); etc. It follows, using ([S]), that the Mattis Hamiltonian in ([T]) can be written 
in terms of free bosons, and this is the key step towards the exact solution of the Mattis 
model 1261. 



The gauge fields are quantized as usual by a partial gauge fixing Aq = 0, postulating 
the canonical commutation relations (see below (|2])), and imposing the Gauss law operator 
constraint on the Hilbert space [33]. However, the corresponding quantum Hamiltonian is 
not obtained as the straightforward quantization of the minimally coupled classical Hamil- 
tonian: due to the anomalous commutators in (jlj), the Gauss law operators G[x] would 
in this case no longer commute with the Hamiltonian, and the theory would thus not be 
gauge invariant. The remedy of this problem is to introduce a manifestly gauge- invariant 
normal-ordering prescription. We use the point-splitting method pioneered by Schwinger 
[32] : start with the gauge-invariant expression 

V^t^(x - ee./2) e'^'^I-i%M^^ie.),,^^^^^^ ^ ^^^/2), (7) 

which includes a line integral of the gauge field. We define our gauge-invariant fermion 
normal ordering of bilinears, liplgip^go, as the limit e — )■ of ([7]) after first subtracting off 
its singular part r/{2nide) (see Equation(4.9) in [26]). The result is 

oipl^ip^/o = J^^, + reoAJ{2-Ka), (8) 

and similarly 

d^x li)l,r {-ids + e^As) i),/o = I d^x (^:^,\,r(-i9,)^,_,: +reoA,J,,, + ^^s) (9) 

(computational details will be provided elsewhere). An important feature of ([9]) is the 
"bare" photon mass term in the second line; as noted, this is a direct consequence of gauge- 
invariant normal-ordering. One can verify that (j8]) and (|9]) are gauge- invariant expressions. 

The Hamiltonian and the Gauss law operators of the gauged model are now bosonized 
using the above results: 



H=]^j d^xl{vFY. [(1 - 7i)(n, - euAsf + (1 + 7i)(5.<l> 

+72(5.<l>.)(9_,<l>_. 



(10) 



,2 d2 \ X 



I- e; + El + c'B 

and G[x\ = J d'^xxJ2s=±'^si—Es + cr^s), with the renormalized charge Cr = eQ/{-\fna). 
Note that the scaling of the model parameters and boson fields are such that our gauge 
theory model remains well-defined in the UV limit a — )■ 0^. The charge renormalization 
Co ~^ ^R shows that photons can have stronger influence on physical properties than what 
superficial arguments might suggest. Under a gauge transformation, Ag ^ Ag + dgXi 
Es — 7- Eg, $s —7- $s and H^ — > H^ + crOsX, such that H^ — crAs, and thus H in flTOj) . are 
gauge invariant. 



The Hamiltonian in ( ITOl) is quadratic in boson operators and can therefore be diag- 
onalized by a Bogoliubov transformation. To this end, we perform a Fourier transfor- 
mation, -E's(x) — 7- Es[p), and define longitudinal- and tranverse fields by |p|£^i(p) = 
ip+-E+(p) + ip_E^{p) and |p|Et(p) = ip+-E_(p) — ip_-E+(p) (similarly for A^, H^, and $s). 
The Gauss law constraint then implies that |p|(— i^^ + cr^l) is zero on the physical space. 
Fixing the Coulomb gauge, A^ = 0, and solving the Gauss law. El = ^r^l, we obtain the 



gauge-fixed Hamiltonian Hgj, describing transverse photons (At) coupled to longitudinal- 
and transverse plasmons ($l and $t, respectively). By straightforward methods, this 
Hamiltonian can be written in the diagonal form -f^g./. = Eq + X],=i 2 3 Xlp '^i(p)^!(p)^j(p) 



with standard boson operators, i.e. [bj{p),bj,{p')] = 5jj'6p^p', etc. (see e.g. Appendix C in 
[26] for further details), and the groundstate energy Eq. The exact dispersion relations 
ujj{p) are computed from the eigenvalues of a certain 3x3 matrix. We obtain the following 
characteristic polynomial of this matrix whose zeros A = Xj are equal to Ujlp)"^ 

A(A - e^ - c2|pn(A - e^ - vl\p\^) + IpI^SV (4(A - c^IpP) - c202) (11) 

with 4 = vUl - 7i)(l - 71 ± 72)/2, v^ = vUl - 7?), e = v/^^(l-7i)|e^J, and S = 
\sm{2ip)\ = |2p+p_|/|pp (IpP = p"^ +P-)- We thus obtain two gapped modes with the 
same gap proportional to the renormalized charge: a;i(0) = a;2(0) = 0, whereas the third 
mode is gapless, ^3(0) = 0. Moreover, for c^ vp and f|.|ppS' ^ 1, a;i(p) ~ a/0^ + c^IpP, 



0^2 (p) ~ \/©^"+"yMpP, and 




This shows that ui gives the energy spectrum of the dressed transverse photon, whereas 
U2 and Go's give those of the dressed transverse- and longitudinal plasmons, respectively. 
Remarkably, the behavior of the last mode is qualitatively different for cr = and e/j 7^ 0: 
in the former case, oosip) ~ \p\Sy/v+V-/vF, which is equal to the lowest-energy mode of 
the Mattis model [26], whereas in the latter case, uj^^p) ^ |ppS'cf_/9 for p < Q/c. Thus 
the photons can affect the low-temperature thermodynamical properties of the system, no 
matter how small IcrJ 7^ 0. For example, we found that the temperature (T) dependence of 
the specific heat at low T is linear for e^ = 0: C^ oc T, but there are logarithmic corrections 
for en ^ 0: C^ oc Tln(To/T) with Tq = {n/ayc^/ {1.4230). 



We also studied the magnetic field response to an external current J^ with J° = 
0, i.e., we computed the linear response function K{ijj,p) in the relation {B{u,p)) = 
i^(a;, p)|p| Jt(w, p) + 0{J^) (IpI Jt is the Fourier transform of 9+J_ — 9_J+). The ex- 
act result is 

^ _ a;^(a;^ - 6^ - v].\p\') + vlS'lpWe' + t;i|p^|) 

^'^'^^ (-a;U^i(p)^)(-^U^2(p)^)(-u;U^3(p)^) ^ ^ 

with w^ = (w + iO^)^ and ouj{p) the dispersion relations given above. Expanding K{u,p) 
in partial fractions and inserting the long-distance approximations of u}j{p) given above, 
we obtain for c^ vp, 

{B{U, p)) ^ _^2 ^c^\p\2^Q2 \P\MuJ, p), (14) 

ignoring 0( J|.)-terms. This proves that there is a Meissner effect with a London penetration 
depth Xl = c/6. The question if our model has other features of a superconductor, e.g. a 
fermion condensate, is under investigation. 



The exact solution of the model presented in this paper provides a reliable starting 
point to explore the physics of 2D correlated fermion systems. If the gauge field in this 
model is to be interpreted as a physical electromagnetic field, one should extend it to three 
dimensions. Furthermore, if the model is to describe correlated fermions in real materials, 
one should introduce spin degrees of freedom, and investigate if the gapped fermions affect 
the ones in the nodal regions. Exactly solvable extensions of our model addressing these 
remarks will be presented elsewhere. 
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